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r~| ' Abstract 

C^ I We show that the time-evohition of the wave hmction describing the macroscopic variations 

of the pair density in BCS theory can be approximated, in the dilute hmit, by a time-dependent 
Gross-Pitaevskii equation. 

1 Introduction and main result 

OO , We consider systems of fermions interacting through a two-body potential admitting a bound state 

^ . I with negative energy. At very low temperature and density, the fermions are then expected to form 

^!j ■ tightly bounded pairs, which behave hke bosons and produce a Bose-Einstein condensate. From 

^sg . the mathematical point of view, it would be very interesting to establish the validity of this picture 

starting from first principle many body quantum mechanics. This seems however a very difficult 
task. Recently, this problem has been studied in [7], starting from the Bardeen-Cooper-Schrieffer 
(BCS) approximation of many body quantum mechanics. In this paper, the authors show that at 
zero temperature the macroscopic variations in the pair density are described (to leading order) by 
the Gross-Pitaevskii energy functional, corresponding to a Bose-Einstein condensate of fermion pairs 
with an effective repulsive interaction. In the present paper, we analyze the same regime, but from 
a dynamical point of view. We consider the time evolution, governed by the time-dependent BCS 
equation, of an initial BCS state with sufficiently small energy. The energy condition guarantees that 
fermions form pairs, and it fixes the microscopic structure of the pair density (see Proposition 11.21 
below). We assume that, at time i = 0, the macroscopic variations of the pair density are described 
by a wave function ip € L^(]R^). Then we show that, at time t 7^ 0, on the macroscopic scale, the pair 
density is described by a new wave function ipt € L^(R^), given, to leading order, by the solution of 
the time-dependent Gross-Pitaevskii equation with initial data V'i=o = ^• 

As pointed out above, our analysis is relevant for dilute systems at or close to zero temperature. 
If instead one considers the regime close to the critical temperature, the macroscopic variations of the 
pair density are described by the Ginzburg-Landau energy functional; this has been recently proven. 
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starting again from the BCS approximation, in [3], making use of earlier works [U [5l [TO] where, 
among others, the existence of a unique critical temperature was established for systems interacting 
via a general class of two body interactions. Let us now summarize the physical picture. In the 
seminal paper [12] Leggett suggested that the BCS functional, introduced by Bardeen, Cooper and 
Schrieffer in [T], serves as a valuable description of a large class of Fermi systems, with interactions 
ranging from weak two-particle potentials up to rather strong interactions allowing for bound states. 
For weak attractions, fermions tend to form Cooper pairs, living on a much larger scale compared 
with the mean particle distance, and displaying superfluidity (superconductivity). In the case of 
strong potentials, the fermions are forced into tightly bound pairs behaving as bosons and forming a 
Bose-Einstein condensate. This picture was later extended to positive temperature in [H]. Whereas 
the BCS functional is supposed to remain valid in the whole crossover region from weak to strong 
coupling, two distinct parameter regimes are of particular interest, as pointed out in J13l [3], One is 
the low density and low temperature limit which, for strong coupling, leads to the Gross-Pitaevskii 
theory (the BEC regime). The other one is the limit close to the critical temperature, where the 
macroscopic variations are captured by the Ginzburg-Landau theory. For nice reviews on the physical 
background of this subject we refer to [15^ [2]. Whereas the emergence of these effective models on 
the macroscopic scale was proven in [HI [7] in the static case, the goal of the present paper is to settle 
the evolution problem on the BEC side. A very interesting open problem, which is not addressed in 
the present paper, consists in understanding (starting from the BCS approximation or, even more 
ambitiously, from many-body quantum mechanics) the emergence of a time-dependent Ginzburg- 
Landau type equation for the description of the macroscopic evolution in the regime close to the 
critical temperature. 

In BCS-theory, the state of the fermionic system is described by a self-adjoint operator T : 
L2(r3) e l2(m3) -^ L2(r3) l2(r3), satisfying < P < 1, defined by the 2 x 2 matrix 

Z " ) , (1.1) 

a 1 — 7 y 



where 7, a : L (M ) — )■ L (R ) will be described in terms of their integral kernels 7(2;, y), a{x,y). 
The bar indicates complex conjugation, i.e. a{x,y) = a{x,y). The fact that F is hermitian implies 
that 7 is hermitian and that a is symmetric, i.e. ■j(x,y) = 7(2/, x) and a{x,y) = a{y,x). Moreover, 
the condition < F < 1 implies F^ < F and thus < 7 < 1 and aa < 7(1 — 7). There are no spin 
variables in F. We implicitly use SU{2) invariance of the states. In fact, one has to imagine that 
the full, spin dependent. Cooper pair wave function is the product of a{x,y) with an antisymmetric 
spin singlet. 

We are going to assume particles interact through a two-body potential V admitting a negative 
energy bound state. More precisely, we will assume that V satisfies the following assumptions. 

Condition 1. We assume V G L^(]R3) nL°°(]R3), with V{—x) = V{x). Moreover, we assume that 
— 2A -|- V has an isolated ground state energy —E^ < with a unique non-negative normalized 
ground state oq, satisfying (— 2A -|- V)aQ = —Ebao. Under these assumptions, ao{—x) = ao{x), 
ao G L1(M3) n L°°(M3) and {xl'^aoix) G L" 



We are interested in the ultradilute regime, where A^ ^ 1 fermions move in a volume of order A^^ (we 
consider a system defined on R^; the length-scale is determined here by the choice of the initial data). 
The particles are also subjected to an external potential, which we assume to be weak and slowly 
varying (it varies only over macroscopic scales of order A^). It is useful to rescale lengths, introducing 



macroscopic coordinates. Putting h = 1/N ^ 1 (so that the expected number of particles in the 
system is I//1), the BCS energy functional in macroscopic units is given by 

^BCs(r) = Tr {-h^A + h^W) 7 + ^Jdxdy V (^) Hx,y)\^. 

For a formal derivation of the BCS functional from quantum mechanics, see e. g. [S", Appendix A]. In 
[7], in a slightly different setting (with periodic boundary conditions), it was proven that, for h <^ 1, 

inf ^BCs(r) = -§ + /i „ inf ^ ^GP(V') + o{h) 

0<r<i;Tr7=i/fc 2/1 ||Vll2=i 

with the Gross-Pitaevskii energy functional 

£gpW= [dx Q|VV'(x)|2 + 2M^(x)|V'(x)|2+5|^(x)r 

and the coupling constant 

g=j^ldq\ao{q)\U2q^ + E,). (1.2) 

It was also shown in [7] that the pair density of every approximate minimizer of the BCS functional 
has the form 

oi{x,y) = -^ijj(^—j ao(^—j +^{x,y) (1.3) 

where tp is an approximate minimizer of the GP functional, and where ||^||2 < Ch}-'"^ and ||C|[hi — 
Ch~^''^ as h ^ (for comparison, ||a||2 — /i~^'^ and ||a|l//i — h^^'"^, as /i — )• 0). The interpretation 
of the results of [7] is straightforward: at zero temperatures, particles will minimize their energy 
by forming pairs, each having energy —Eb (the main term in the energy, —Ei,/2h, is exactly the 
energy of the l/2/i = A^/2 pairs). The pairs will then behave like bosons and, in the ultradilute limit, 
they will form a Bose-Einstein condensate, with condensate wave function given by the minimizer of 
the GP-functional (the second term in the energy, given by /iinf||^||2=i £'gp(V'); is the energy of the 
condensate) . 

In the present manuscript we consider the dynamics in the ultradilute limit, as determined by 
the BCS theory. We are going to study the evolution of BCS states with energy bounded above by 

SBCsir) <-^ + Ch (1.4) 

in the limit of small h <^ 1. This assumption on the energy implies that the fermions will still 
form pairs or, equivalently, that the microscopic structure of a{x,y) is still given by ao((x — y)/h), 
as in ()1.3p . In other words, (jl.4p guarantees that only macroscopic excitations of the pair density 
are allowed, because microscopic excitations would cost too much energy. Since the relative energy 
AE/E of the excitations we are considering is of the order /i^ , we have to consider times of the order 
l//i^ in order to observe a non-trivial dynamics. With respect to rescaled time, the (Hamiltonian) 
time evolution generated by the BCS functional <?bcs is governed by the equation 



I 



h^dtTt = [Hr„rt] (1.5) 



for the time-dependent operator Tt : L'^iM.^) ^^(M^) -> L'^(M.^) L'^{R^). The self-consistent 
Hamiltonian ffrt is given here by 



^r. 



-h'^A + h^W Vat 

Vat /i^A - h^W 



A similar equation emerging in the study of the dynamics of stars was recently studied in [9]. As 
discussed in [9] the proof of existence of local in time solutions of (jl.Sp is standard. Since local 
solutions preserve the energy, and since the potential is assumed to be bounded, one immediately 
obtains global well-posedness for (jl.Sp . for initial data in the energy space. Eq. (jl.Sp translates into 
the following two coupled nonlinear equations for the kernels 7((x,y) and at{x,y): 



ih^dt at{x, y) = (-h^A, + h^W{x) - h^Ay + h^W{y) + V ( ^^ ) ) at{x, y) 



at 

' X — y 

—jj^ty-^y' (1.6) 



z — V \ I / z — X 

dz-ft{x, z)V ( — — j at{z, y) - dz-ft{y, z)V ( —^ ] at(x, z) 



with the operator Gat defined by its integral kernel 



Gat (x, y) =i I dza{x, z)a{y, z) (V ( ——- ] -V 



h J \ h 

Note that Tr7t is preserved, because Tr Gq,j = for all t G M. In the following we will focus on the 
equation for at- It is useful to represent the kernel at using center of mass and relative coordinates. 
For this reason, we define 

at{X,r):=at{X + r/2,X-r/2) 

which means that at{x, y) = at{{x + y)/2, x — y). From (jl.Sp . we obtain then the following equation 
for the evolution of at'- 

ih^dtatiX, r) = f-y Ax - 2/1^ A, + V{r/h)] at{X, r) 

+ h^ (WiX + r/2) + WiX - r/2)) 5t(X, r) 

(1.7) 
dz-it{X + r/2,X + z - r/2)V (z/h) 5t(X - r/2 + z/2, z) 



dzjt{X - r/2, X -z + r/2)V [z/h) at{X + r/2 - z/2, z) 

We rewrite the equation (|1.7p in integral (or Duhamel) form: 

1 /•* 

at{X, r) = U{t)ao{X, r) + - / ds U{t - s) {W{X + r/2) + W{X - r/2)) 5,(X, r) 

1- Jo 

-^ [ dsU{t-s) f dz^s{X + r/2,X + z-r/2)Viz/h)asiX -r/2 + z/2,z) (1.8) 

-— 2 / dsU{t-s) dz-fs{X -r/2,X - z + r/2)V{z/h)as{X -r/2- z/2,z) 
where we defined the evolution 

with Ux (t) = e*(*' ^-^ -^ denoting the free evolution in the center of mass coordinate X and 

^^(^) ^ ^-it{-2Ar+il/h^)Vir/h)) 

the evolution in the relative coordinate r. Observe that 

Ur{t)ao{r/h) = e''^''/''\o{r/h) (1.10) 

We are now ready to state our main results. 



Theorem 1.1. Suppose Condition 1 above is satisfied, and assume W € -ff^(M^)nL°°(M^). Consider 
an initial BCS state < Fq < 1 with Tr'jQ < C/h, and 

£BCs{To)<-^Tr^^ + Ch (1.11) 

Let at he the solution of the integral equation \1.^) . and set, for arbitrary t € M, 



MX) -^J^ I drao{r/h)at{X,r) (1.12) 

Then we have ||V't||//i < C, uniformly in h and in t € M, and 

-itE^/h^ 

at{X,r) = ^^ MX)ao{r/h) + e"**^''/'^ UX,r) 

with ||Ct||^2™3xM3) ^ dh and ||'?t||^i™3xM3) ^ Ch~^. Moreover, if Lpt{X) denotes the solution of the 
Gross- Pitaevskii equation 

idm = -^Aift + 2Wipt + 2g\ipt\^ipt (1.13) 

with the initial data ipt=Q{X) = 'iPq{X) and with the coupling constant g as in il.2\) . we have 

llV'i- 9^*112 <C/ii/Vl*l (1.14) 

for constants C,c > 0. 
Remarks: 

• The theorem impHes that, in good approximation, the pair density at is given by the product 
of the ground state wave function oq (varying on the microscopic scales, describing the internal 
degrees of freedom of the fermion pairs and remaining constant in time) and the solution of the 
Gross-Pitaevskii equation, varying on the macroscopic scale and describing the distributions of 
the pairs in space. 

• We only assume that the number of particles Tr7o is bounded above by C/h. The interesting 
regime is the one with Tr7o of the order 1/h. It is easy to see that HV'ilb — hTr^o; hence, if 
Tr7o <^ l//i, the convergence (J1.14p is less interesting (it becomes trivial if Tr7o < l/h^'"^). 

• The Gross-Pitaevskii equation ()1.13p is a defocusing non-linear Schrodinger equation. It is 
known to be globally well posed in the space i:f^(M'^). In other words, for any initial data 
(po € -ff^(M^), there is a unique solution ift G C{{—oo,oo),H^{M.^)). Since the nonlinearity is 
defocusing and the external potential is bounded, it is easy to see that the H^ norm of ft is 
uniformly bounded in time; i.e., there exists a constant C, depending only on the iJ^-norm of 
the initial data, with llv^tH//! < C for all t € M. 

• The fact that tpt represents a pair of fermions is nicely reflected in the Gross-Pitaevskii equation. 
The factor 1/2 in front of the kinetic term corresponds to the total mass m = 1 of two fermions 
(we use units where each fermion has mass 1/2). The factor 2 in front of the external potential 
W corresponds to the total charge e = 2 of the two fermions. 



• Our analysis can also be extended to particles in small external magnetic fields, varying on 
the macroscopic scale. The magnetic field acts on the particle through a minimally coupled 
magnetic potential, which, in macroscopic units, has the form hA{x) = h{Ai{x) , A2{x) , A^{x)) . 
The limiting Gross-Pitaevskii equation ()1.13p is replaced in this case by 

idm = ^{-iV - Af^t + 2W^t + 2g\iptf^t (1.15) 

The proof of Theorem 11.11 extends to a derivation of (jl.lSp , under the assumptions that A is 
continuous and decays at infinity and that the spectrum of (— iV — A)'^ is absolutely continuous 
(these conditions guarantee the validity of Stricharzt estimates similar to those discussed in 
Appendix |A] for the evolution generated by the magnetic Laplacian; see [6]). 

A crucial ingredient of our analysis is the following energy bound established in [H [7], which we 
adapt here to our setting. 

Proposition 1.2. Suppose Condition 1 is satisfied, and W € L°°(K'^), and consider a BCS state 
< r < 1 with Tr-f < C/h and 

SBCs{T)<-^Tr^ + Ch. 

Then Try"^ < Tr^j — iya) < Ch for C > independent of h (from 7(1 — 7) > aa we find that 
7 — aa > 7^ and in particular that 7 — aa is a non-negative operator). Moreover, with 

^(X) ■■=\ I drao (r/h) a{X,r) (1.16) 

we have HV'IIhi ^ C, uniformly in h > 0, and 

a{X,r) = ^i;{X)ao{r/h)+^{X,r) 

with the bound ||^||^2 ^ Ch and ||^||^i < Ch^^. In particular, since the number of particles Trjt 
and the BCS-energy Sscsi^t) o,fe preserved by the BCS time-evolution, the solution Tt of the time 
dependent BCS equation with an initial BCS state < Fq < 1 with TrjQ < C/h and 

£bcs{^q) < — ^^^^70 + Ch 
is such that Tr^f < Tr{jt — «*«*) ^ Ch and such that, with the definition 

V'i(X) := / drao{r/h)at{X,r) , 

one has WiptWH^ ^ C, uniformly in h and in t, and 

-itEt/h^ 

atiX,r) = ^^ MX)ao{r/h) + e-''^"/'' Ct{X,r) 

with the error ^t satisfying ||^t||^2 ^ Ch and ||Ct||^i < Ch~^ . 

Remark: We will also use the notation ^s{x,y) = ^^((a; + y)/2,x — y) (so that as{x,y) = 
{l/h'^)ips{{x + y)/2)ao((x — y)/h) + ^s{x, y)). The bounds for ^s immediately imply similar bounds 
for is- 



Proof. We start by noticing that 

Eh 
Ch > ^BCs(r) + ^Tr7 



= TT(-h^A + h^W + ^)j + ^jdxdyV{{x-y)/h)\a{x,y)\^ 
> Tr (-h^A + y) 7 + ^ f dxdyV{{x - y)/h)\a{x,y)\^ - h^\\W\\o^TTj 
Since Tr7 < C/h, and since 7 — aa > 7^ > 0, we find 

Ch > Tr (-h^A + Y ) (^ ~ "") + ^^ ( -^^^ ^2)^^^\ I "^^^^ ^^^"^ ~ ^)/^) l"(^' y)l^ 
> C Tr (7 - aa) + Tr ( -/i^A + -^ j aa + - / dxdy y((x - y)/h)\a{x, y) ^ 

Next we introduce center of mass coordinate X = (x + y)/2 and relative coordinate r = x — y. Since 

Ar, + Ay = (l/2)Ax + 2Ar, we conclude that 

TrAaa = / dxdy A^a{x,y) a{x,y) = ^ / dxdy {A^ + Ay)a{x,y)a{x,y) 
= I dXdr (jAx + A,, j S{X, r) 5(X, r) 



Therefore 

h^ 



Ch > CTr (7 -aa) + — / dXdr |Vx5(^,r)p 

+ jdX (a{X, .), ("-/i^A, + W{r/h) + ^^'j 5(X, 



(1.17) 



This implies first of all the bound Tr 7^ < Tr(7 — aa) < Ch (the first inequality is a general property 
of BCS states, since 7(1 — 7) > aa). Now we set 

a{X, r) = ^i;{X)ao{r/h) + ^{X, r) (1.18) 

with ijj as defined in (fTTB]) . From (fTTTl) (using that (-A + (l/2)y)ao = {Eb/2)ao), we find 

/i2 



Ch > C7Tr (7 - aa) + — / dXdr \Vxa{X,r)\^ 

+ j dX (^{X, .), (^-h^Ar + \vir/h) + ^) CiX, .) 



(1.19) 



Eq. (jl.lSp implies that, for (almost) every X € M^, S,{X,.) is orthogonal to oq. Hence, if k > 
denotes the spectral gap between the ground state of —A + (l/2)y and its first excited state, we 
deduce from (J1.19p that 

j dXdr\l{X,r)\^ < Cn'^h 



The L^-norm of ip can be estimated using the very definition (|1.16p . By Cauchy-Schwarz, we find 

||V'll2 - TI / dXdridr2ao{ri/h)ao{r2/h)a{X,ri)a{X,r2) 
< ^ f dXdridr2\aoir2/h)\'^\a{X,n)\'^ < /i||5||i < C 

because ||a||2 = Traa < Tr7 < C/h. To bound the L^-norm of VV' we use again Cauchy-Schwarz 
to estimate 

llVV'lli < ^ f dXdndr2aoin/h)aoir2/h)Vxa{X, n) Vxa{X, ra) 
< -^ f dXdridr2\ao{r2/h)\'^\VxaiX,ri)\'^ < hWVxag < C 
where we used (J1.19p to bound |lVx5||2 < C/h. Since 



^V'(X)ao(r//i) 



' = -1 I dXdr\V^{X)\^al{T/h) = ""°"^j|^^"^ < C/i-i 



the bound HVxaJli ^ C/h also implies that HVx^lJ^ < Ch ^ . On the other hand, since V ^ L^ 
we have -2h^Ar + V{r/h) > -2h^Ar - C and therefore, by (IlTqI) . 

Ch> I dX (I{X, .),[-2h^ Ar - C]l{X, .)\ =2h^ j dXdr\Vrl{X,r)\^ -C\\i\\l 
Since H^lH < Ch, we immediately conclude that HVrClll — Ch^^. D 

2 Proof of Theorem 11.11 

We start from the definition (J1.12p for the wave function ^pt, that is 

i'tiX) := 7 / draoir/h)atiX,r), 

and we insert the integral equation (jl.Sp for at- We find 



MX) 

" h 

^itEb/h- 



drao{r/h)U{t)ao{X,r) 



+ drao{r/h) / dsU{t - s){W{X + r/2) + W{X - r/2))asiX,r) 

ih J Jo 

itEb/h? r r^ r 

drao{r/h) / ds U{t - s) dz 7,(X + r/2, X + z- r/2)V{z/h) 



I 



h^ 







x5,(X-r/2 + z/2,z) 



itEi,/h'^ p ft 



drao{r/h) / ds U{t - s) / dz-fs{X - r/2, X -z + r/2)V{z/h) 



ih^ J Jo 

xa,{X + r/2-z/2,z) 



We use now the decomposition (jl.Op of the evolution U{t), saying that U{t) = Ux{t)Ur{t). Since 
we integrate over r we can let Ur{t) act on the left in the form of its adjoint and use (jl.lOp . more 
precisely U*{t — s)aQ = ^-'^{i-s)^b/h q,^^ such that the phase in the t- variable cancels. We conclude 
that 

MX) 

= Uxit)MX) + ^ I dsUxit-s) fdrao{r/h){W{X + r/2) + W{X-r/2))e''^>'/'''s,{X,r) 

-T^ f dsUxit-s) f drdzaoir/h)-fsiX + r/2,X + z-r/2)Viz/h) 

xe''^''/^^as{X + {z-r)/2,z) 

- — 3 / dsUx{t-s) I drdzao{r/h)-fs{X -r/2,X -z + r/2)V{z/h) 

xe''^»/''''as{X + {r-z)/2,z) 

(2.20) 

Notice that we can see already the Gross-Pitaevskii equation appear. To leading order the integration 
over the function ao{r/h) forces r to be close to in the above integrals, such that one can see the 
term with the potential 2W{X) emerge. For the non-linear part we will simply use the fact that 7^ 
is essentially a^Us such that the last term has to include IV'PV' to leading order. 

By Proposition 11.21 the second term on the r.h.s. of the last equation ()2.20p can be decomposed 
as follows: 



J- f dr ao{r/h){WiX + r/2) + WiX - r l2)y'^^l^''a,{X, r) 



"^'^^^ ' dral{r/h){W{X + r/2) + W{X - r/2)) 



+ ^ f draoir/h){WiX + r/2) + W{X - r/2))^,(X, s) 

2WiX)MX) 
+ MX) [ dral{r)iWiX + hr/2)-W{X)) 



+ MX) / dra'o{r)iWiX - hr/2) - W{X)) 
+ ^ f drao{r/h){WiX + r/2) + WiX - r/2))^,(X, r) 
=: 2W{X)MX) + Gf}{X) + Gg(X) + Gf}{X) 
where we used the normalization Ijaolb = 1- We observe that 



\G^'>{X)\<{h/2)\i^,{X)\ dratir 



/ dKr-VW{X + hKr/2) 
Jo 



< ih/2)\ij,{X)\ I dr I dKaQ{r)'^\r\\VW{X + hKr/2)\ 



Therefore 



I,s 116/5 



< 



C7/i6/5 f I dr\r\al{r) 



6/5 

dr I dK\r\a^{r)\VW{X + hKr/2)\ 



1/5 i- r-l 

dXdr / dK\r\al{r)\iPs{X)f/^\VW{X + hKr/2)f/^ 
Jo 



where we used Holder inequahty with exponents 6 and 6/5 (with respect to the measure dK(ir|r|a^(r)). 
From the assumptions on W, Condition 1 and Proposition !!. 21 (which imphes that 11-0^11^1^1 is uniformly 
bounded in h and in s), we find 



GWl<Ch^/^ 



1 6/5 



dXdr\r\alir)\i;s{X)f 



2/5 



dXdr\r\alir)\VW{X)\'' 



3/5 



< ch^/^UsWHiW^wwl^^ < ch^/^ 



The term G2 l can be bounded analogously. To bound G^l, on the other hand, we observe that 



ig: 



(0) 
3,s 


,6/5 
l6/5 


< 


1 

/16/5 






< 

< 


1 

/16/5 

c 



drao{r/h){W{X + r/2) + W{X - r/2)%{X, r) 

1/5 



/l3/5 



dX 
drao{r/h) 
dXdrao{r/h)\W(X + r/2) + Pr(X - r/2)f 

dXdrao{r/h)\l{X,r)\^ 



6/5 



dX(irao(r//i)|VF(X + r/2) + Ty(X - r/2)f/^\Cs{X,r)f/^ 

2/5 



3/5 



< 



C/l^^^||W^|l3^^||6ll2^^<C/l6/^ 



for a constant C depending on ||ao||i) on ||"o||oo! and on ||l^||j|^i. Here we used again Proposition 
11.21 to bound ||^s||2 < Gh}'"^ and the assumptions on W. We conclude that 



3 
i fdrao{r/h){WiX + r/2) + WiX - r/2))e^^^''/'^'5,(X,r) = 2W{X)MX) + ^GfliX) 
'' i=i 



(2.21) 



.(0) 



where ||Gj jUg/s < Gh, for all j = 1,2,3 and for a constant C > independent of h and of time s. 
We consider now the third term on the r.h.s. of ()2.20p . Using Proposition II. 2^ we find 

A:= - i /" drdzaQ{r/h)-is{X + r/2, X + z- r /2)V{z/h)e''^''/^^as{X + {z - r)/2, z) 



1 
1 



drdzaQ{r/h)-is{X + r/2, X + z - r/2)(yao)(z//i)V's(^ + {z - r)/2) 



^ / drdz ao(r//i)7,(X + r/2, X + z - r/2)F(z//i)e^^-^''/^ ^.(X + (z - r)/2, z) 

Next, we write 7^ = a^as + (7s — SsCts). Here ^s — oisOis is a positive operator (since 7^(1 — 7^) > agOts 
implies that 7^ — agCHs > 7^), with Tr (7^ — asOtg) < Gh (compared with Tra^as ~ Gh~^); we will 
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use this fact to show, in Lemma 13.21 that the contribution from (7^ — Os^s) is small, in the limit 
h^O. We find 

A= --^ drdzdw Qo(r//i)a^(X + r/2, X + r/2 + w)asiX + r/2 + w,X + z- r/2){Vao){z/h) 

xi;s{X + {z-r)/2) 
_ J_ /" drdz ao(r//i)(7. - asas)iX + r/2, X + z- r /2)V{z/h)^,{X + {z - r)/2)ao{z/h) 

-^ f drdz ao{r/h)js{X + r/2, X + z- r /2)V{z/h)e''^''l^%{X + {z - r)/2, z) 

We express again the kernels a^ using center of mass and relative coordinates. We find 
1 



A= - — I drdzdw ao{r/h)as{X + (r + w)/2, w)as{X + {z + w)/2, r + w- z){Vao){z/h) 

X Vs(X + (z-r)/2) 
1 



^5 drdz ao{r/h){-fs - asas){X + r/2, X + z- r/2)V{z/h)MX + {z - r)/2)ao{z/h) 

-^ I drdz ao(r//i)7,(X + r/2, X + z- r/2)V{z/h)e'''^^/^^UX + {z - r)/2, z) 

Finally, we use again Proposition 11.21 to rewrite the kernels a^. We find that 

A= --^ drdzdw ao{r/h)e''^''/''''lp,{X + (r + z)/2)ao{w/h)as{X + {z + w)/2, r + w - z) 

x{Vao){z/h)MX + {z-r)/2) 



. drdzdw ao{r /h)e''^>'/'' ^s{X + {r + w)/2,w)as{X + {z + w)/2,r + w- z) 

x{Vao){z/h)MX + {z-r)/2) 
^ f drdzao{r/h){js - asas){X + r/2, X + z- r/2)(yao)(^/MV's(X + {z - r)/2) 

^ / drdzao{r/h)-f,{X + r/2, X + z- r/2)V{z/h)l{X + (z - r)/2, z) 
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and we conclude 
1 



A = - — / drdzdw ao{r/h)ao{w/h)i;^{X + (r + z)/2)i;s{X + {z + w)/2)ao{{r + w- z)/h) 



1 



>i{VaQ){z/h)^s{X + {z-r)/2) 
drdzdw ao{r/hytp^{X + (r + z)/2)ao{w/h)Cs{X + {z + w)/2, r + w - z) 

x{Vao){z/h)MX + {z-r)/2) 



-y / drdzdw ao{r/h)^s{X + (r + w)/2, w)aQ{{r + w- z)/h)^siX + (z + w)/2) 

x(yao)(z//i)V.(X + (z-r)/2) 

-T I drdzdw ao{r/h)is{X + (r + w)/2, w)ls{X + {z + w)/2, r + w - z) 

x{Va^){z/h)^s{X + {z-r)/2) 
drdzao{r/h){js - asas){X + r/2, X + z- T/2){Vao){z/h)i^s{X + {z - r)/2) 



1 
1 



drdz ao{r/h)js{X + r/2, X + z- r /2)V{z/h)is{X + {z - r)/2, z) 

5 



-.M{X) + Y,Gfl{X) 



(2.22) 



We decompose further the main term M(X), writing 

M{X)= — Tg drdzdw aQ{r / h)ao{w / h)ao{{r + w — z)/h){Vao){z/h) 

X i^,{X + r/2 + z/2)i;s{X + z/2 + w/2)i;s(.X - r/2 + z/2) 
drdzdw aQ{r)ao{w)ao{r + w — z)(yao){z) 

X V^,(X + /i(r + 2)/2)V's(X + h{z + u;)/2)Vs(^ + Hz - r)/2) 
4's{X) I drdzdw aQ{r)aQ{w)aQ{r + w — z){VaQ){z) 
X i^,{X + /i(r + z)/2)i)s{X + h{z + w)/2) 
drdzdw ao{r)ao{w)ao{r + w — z)(yao){z) 

X V^,(X + hir + z)/2)iJsiX + /i(z + u;)/2) {iPsiX + /i(z - r)/2) - Vs(^)) 
'(/'s(-^)^ / drdzdw ao(r)ao(tf)ao(^ + "W^ ~ -2)(^"o)(-2)^s(-'^ + /i(?' + -z)/2) 

i^s{X) / drdzdtt7Qo(r)ao(w')ao(r + tf — z)(yQo)(z) 

X V^,(X + hir + z)/2) (V-slX + h{z + ti;)/2) - MX)) 
/ drdzdw Qo(r)ao(tf)ao('^ + tf — z)(yQo)(z) 

X i^,{X + /i(r + z)/2)iJs{X + /i(z + w)/2) (^,(X + /i(z - r)/2) - ^,(X)) 
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In the first term, we isolate again the main contribution. We find 

M{X)= - \'iljs{X)\^iJs{X) drdzdwao{r)ao{w)ao{r + w-z){Vao){z) 

-i{js{Xf / drdzdw ao{r)ao{w)aoir + w - z){Vao)iz) (^^{X + h{r + z)/2) -^^(X)) 

— ips{X) I drdzdw aQ{r)aQ{w)ao{r + w — z){VaQ){z) 

X ^,{X + h{r + z)/2) {ijs{X + h{z + w)/2) - ^s{X)) 

— / drdzdw aQ{r)aQ{w)a{){r + w — z){VaQ){z) 

X V^,(X + h{r + z)/2)i)s{X + h{z + w)/2) (^,(X + h{z - r)/2) - V',(X)) 
We observe that 



— / drdzdw aQ{r)aQ{w)aQ{r + w — z){V aQ){z) 

= - dz{ao *ao*ao) {z){Vao)iz) = ——3 / dp\ao{p)\'^{2p'^ + E^) = g 

where we used the eigenvalue equation (— 2A + V)aQ = —EhOo to write Voq = —Ef,ao + 2Aao and 
therefore Vao{p) = —{Eb + 2p'^)ao{p). We conclude that 

3 

.(2), 



with the three error terms 



M{X) = g\^j{X)Mx) + Y^G^ll{X) 



G['l{X) = -iPsiXf / drdzdw ao{r)ao{w)ao{r + w- z){Vao)(z) (V,(X + h{r + z)/2) - V',(X)) 

G2 s{X) = —tpsiX) / drdzdw aQ{r)ao{w)aQ{r + w — z)(ya()){z) 

X ^,(X + h{r + z)/2) iiJsiX + Kz + w)/2) - V^,(X)) 

G^l{X) = — I drdzdw aQ{r)aQ{w)aQ{r + w — z){VaQ){z) 

X V^,(X + /i(r + z)/2)Vs(^ + h{z + 'u;)/2) (Vs(^ + h{z - r)/2) - ^s{X)) 

(2.23) 

Together with (I2.22|) . this implies that 

A = --^ / drdzaQ{r/h)-is{X + r/2, X + z- r/2)V{z/h)e''^^/^^as{X - r/2 + z/2, z) 

= g \ux)?ux) + i2 G^U^) + E ^S(^) 

In Lemmas 13.11 and 13.21 below, we show that HGj-jUg/s < Ch for j = 1,2, ||Gj.;||6/5 < Ch}'"^ for 
j = 3,5, ||G4g||g/5 < Gh^/'^, and ||Gj ,, II2 < Ch for all j = 1,2,3, for a constant C independent of h 
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and of s. In particular, HG^'^He/s < Ch^'"^ for all i = 1,2 and j's. Similarly, one can write the last 
term on the r.h.s. of (12.20p as 

-p / drdzaQ{r/h)-is{X - r/2, X-z + r/2)V{z/h)e''^'^'^'as{X + r/2 - z/2, z) 
= g\MX)\'MX) + J2 GfliX) + E ^I'liX) 

where HGj^^Hg/s < Ch}'^ for i = 3,4 and all j's, and for a constant C independent of h and of s. 
Therefore ipj^ and (p:^ imply that 

MX) = Ux{t)MX) + - I dsUx{t-s)W{X)ijs{X) 

1- Jo 

+ ^ I dsUxit-s)\MX)\^MX) + - I dsUx{t-s)Fis) 
Wo Wo 

with the error term F = E,Li GfJ + J2%i{G^ll + Gf}) + E%i{GfJ +Gg) is such that ||F(s)||6/5 < 
Ch}''^ . We compare iptiX) with the solution of the Gross-Pitaevskii equation (J1.13p . which can be 
rewritten in integral form as 

n{X)= Ux{t)MX) + - f dsUx{t-s)W{X)ips{X) + ^ [ dsUx{t-s)\^s{X)\^MX) 

i Jo Wo 

The difference between ipt and ipt is therefore given by 

'iljt-ipt = Ux{t){iJo - ipo) -2i ds Ux{t - s)W{ips - <fs) 

Jo 

-2ig I dsUx{t - s) {liPsl'^ilJs - \fs\^'Ps) -i I dsUx{t-s)F{s) 
Jo Jo 

We use Strichartz estimates (see Appendix [A]) to bound this difference. We find (since we will later 
iterate the next bound, we keep track of the difference of the initial data; at the end we'll set ipo = fo) 



■'t-ft\\L°°{[0,T];L2{ 
< UO - fob + ^Wii^s - '/'s)|Il2([o_T],L6/5(R3)) + 2g \\\i^s\^1ps - I</'s|Vs||l2([o,T],L8/5(M3)) 

+ ll^llL2([0,T],Le/5(R3)) 

< llv5o-^o||2 + 2( / ds\\Wii;s-Vs)\\l/r,] +(/ dslllV'.pV^, -|(/p,|V, "^ 

1/2 



6/5 1 T I / ""'IM'/^sl iys-\-rs\ V^s|l6/5 



I ds\\F{s)\\l/,^ 



< II^^O - V'o|l2 + 2T^/^||W^||3||'0s -(^s||ioo([o,T];L2(M3)) +T^^^ SUp || l^-sl Vs " I'/'s | Vs lle/S 

se[o,T] 
Next, we estimate the difference IUV'sPV's ~ IVsPv^sHe/s in terms of Wips — V^slb, using the a-priori 
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bounds for H^^Hji^i (given by Proposition ll.2p and ||(/?s||j|/i (see the remark after Theorem II. ip . 

IIV't-¥'t||L°°([0,T];L2{R3)) 

< 119=0 - V-Olb + 2T'/^\\W\\3\\iJs - V5s||loo([o,t];L2{R3)) + CT^'^h^/^ 

+ T^/^ sup [|l|'(/'sp(V's -Vs)\\&/b + Wi's'Psi'^s -955)116/5 + WWs^ii's -955)116/5] 

sg[0,T] 

< \W0 - M\2 + CT^'^h^'^ + 2T^I^WhUs - '/5s||l-{[0,T];L2( 

+ rl/2 sup {Us\\l + \Ws\\l) Wi^s-^sh 
se[o,T] 



We choose Tq > such that the factor CT^'"^ appearing in front of the last term on the r.h.s. of the 
last equation is less than, say, 1/2 (observe that Tq depends only on the bounds for the H^ norm of 
ips and of \\W\\ and is therefore independent of h and s). Then we have 

ll^s - V's||l°o([o,t],l2) < 2||^o - V'olb + Ch^/"^ 

for all < T < Tq. For t € [uTq, (n + l)To], we find analogously (since all bounds are uniform in 
time) 

< 8||V'(n-2)To - 95(n-2)To h + ^Ch^'^ + 2Ch^/^ + Ch^'^ < 2^+^^^ - (^olb + 2^+^Ch^/^ 

For initial data ^o = V'O) we conclude that 

for all t G M. This concludes the proof of Theorem ll.li 

3 Control of error terms 

In this section, we estimate the size of the error terms from (j2.22p and (|2.23p . 
Lemma 3.1. Consider the error terms 



Gi'i(X) = -i^siXy / drdzdw ao{r)ao{w)aQ{r + w- z){Vao){z) [i^,{X + h{r + z)/2) - V',(X)) 

G2l{X) = —ijJs{X) I drdzdw a{){r)aQ{w)aQ{r + w — z){VaQ){z) 

X ^,(X + h{r + z)/2) (MX + Kz + w)/2) - ^X)) 

G^i.{X) = — / drdzdw ao{r)ao{w)ao{r + w — z){Vao){z) 

X Ip^iX + h{r + z)/2)^siX + hiz + w)/2) {^X + Hz - r)/2) - ^s{X)) 

(3.24) 

as defined in Ii2.23\) . Then ||G'oj||6/5 < Ch, for a constant C independent of h and of the time s. 
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Proof. With 



iPs{X + h{r + z)/2) - Vs(^) = / dK —MX + hK{r + z)/2) 

Jo dK 

= h dKV^psiX + h{r + z)/2)-{r + z) 
Jo 



we find 



\G\ !.{X)\ < h\ips{X)\ / drdzdw / dK\r + z\ao{r)ao{w)ao{r + w — z) 

X |(yao)(z)| |V^,(X + K/i(r + z)/2)| 

Hence, using Holder inequality with exponents 6 and 6/5 (with respect to the integration measure 

drdzd'wdKao{r)ao{'w)ao{r + w — z)\(yao){z)\), we find 



lof^fj^l < h^/^ ( f drdwdz f dK\r + zfao{r)ao{w)ao{r + w-z)\{Vao){z)\ 
X / dXdrdzdw I dKao{r)ao{w)ao{r + w — z)\[Vao){z)\ 



1/5 



< Ch^l^ 



dXdrdzdw i dKao{r)ao{w)ao{r + w — z)\{Vao){z)\\'il>s{X)\ 
Jo 

dXdrdzdw / (iKao(r)ao(tf)ao(r' + tf — z)|(yQo)(z)| 

JO 



2/5 



X |V^s(^ + /^/i(r + z)/2)|^ 



3/5 



<C/.^/'llV'sf/f 



for a constant C depending on ||ao||i, ||ao||oo, ||?'^ao||i, and ||yao||i (in the second integral, we shifted 
the variable X ^ X — Kh{r + z)/"!). Similarly, we bound 

\G2g{X)\ <h\'ilJs{X)\ I drdzdw I dK\w + z\ao{r)ao{w)ao{r + w — z)\{Vao){z)\ 

X \iPs{X + h{r + z)/2)\ \Vi^s{X + Kh{w + z)/2)\ 
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which gives 



IG^ifJ/l < h^^^ ( [ drdzdw f dK\w + z\^ao{r)ao{w)ao{r + w - z)\{Vao){z) 



1/5 



X / dXdrdzdw / dKao{r)ao{w)aQ{r + w — z)\{Vao){z)\ 

X \i^siX)f^''\MX + hir + z)/2)|6/5 |VV.(X + /./i(w; + z)/2)f/' 



< Ch^/^ 



r-l ll/5 

dXdrdzdw l dKaQ{r)aQ{w)aQ{r + w — z)\[VaQ){z)\\'ilJs{X)\ 
Jo 

dXdrdzdw / dK.ao{r)ao{w)ao{r + w — z)\{Vao){z)\\'ipsiX + h{r + z)/2)f 
Jo 

dXdrdzdw / dKaQ{r)aQ{w)ao{r + w — z)\{VaQ){z)\ 
Jo 

2l3/5 



1/5 



X iVil^siX + Kh{w + z)/2)\' 



<ch'/'mf/^ 



.(2) 



Finally, we bound the term G\ i.. To this end we proceed analogously, noticing that 



\G^l{X)\<h / drdzdw / dK\z — r\aQ{r)ao{w)ao{r + w — z)\{Vao){z)\ 

X |V^,(X + h{z + w)/2)\ iV-sl^ + Hr + z)/2)\ |VVs(^ + Kh{z - r)/2)| 
and therefore concluding that 



\G 



(2)||6/5 
3,s 116/5 



<h' { I drdzdw I dK\z — r\ aQ{r)aQ{w)aQ{r + w — z)\{V aQ){z)\ 
X / dXdrdzdw I dKaQ{r)aQ{w)aQ{r + w — z)\{V aQ){z)\ 



1/5 



X \^s{X + h{z + w)/2)f'' \iJs{X + h{r + z)/2)|6/5 |VV.(X + Kh{z - r)/2)f'' 



< Ch^/^ 



dXdrdzdw j dKaQ{r)ao{w)ao{r + w — z)\{Vao){z)\\ilJs{X + h{z + w)/2)\ 





1/5 



1/5 



dXdrdzdw / dKaQ{r)ao{w)aQ{r + w — z)\{VaQ){z)\\ilJs{X + h{r + z)/2)\ 
Jo 

/•I . 

dXdrdzdw / (i«;Qo(r)ao(^)oo(^ + W7 — z)|(yao)(z)| |V^s(X + k/i(z — r)/2)|^ 
Jo 



3/5 



<C/^^/^IIV'.|l5,/f 



D 
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Lemma 3.2. Consider, from i2.22\) . the error terms 



Gg(X)= -j^ I drdzdwao{r/h)i:,{X + {r + z)/2)ao{w/h)UX + {z + w)/2,r + w-z) 



xiVao)iz/h)MX + iz-r)/2) 



^(1)^ 



G^iiX) 



.{!)/ 



/l5 



G^;^(X)= --y / drdzdwao{r/h)Cs{X + {r + w)/2,w)ao{{r + w-z)/h)'ilJs{X + {z + w)/2) 



x{Vao)iz/h)MX + iz-r)/2) 



drdzdw ao{r/h)^s{X + (r + w)/2, w)^s{X + {z + w)/2, r + w - z) 

x{Vao){z/h)MX + {z-r)/2) 



G\\iiX)= --^ I drdzao{r/h){Vao)iz/h){js - asas){X + r/2,X + z - r/2)MX + (z - r)/2) 



gS(^) 



/l3 



drdzaoir/h)-fsiX + r/2,X + z - r/2)V{z/h)Cs{X + {z - r)/2, z) 



Then \\G, j Hg/s < C/i /or j = 1,2 anc? ||G^- j ||e/5 < Ch^''^ forj = 3, 4, 5, for a constant C independent 
of h and of the time s . 

Proof. We start with the error term G\ !,. We have, using Holder inequahty with the dual exponents 
6 and 6/5, 



\G 



1 



(l)||6/5_ 

l,s 116/5 /),42/5 



dX 



drdzdw aoi^r / h)ao{w / h)(yao){z / h) 



xV',(X + (r + z)/2) i;{X + {z - r)/2)i,{X + {z + w)/2, r + w - z) 



6/5 



< 



1 



/l42/5 



drdwdz ao{r/h)ao{w/h)(yao){z/h) 



1/5 



< 



X / dX drdzdw ao{r/h)aQ{w/h){VaQ){z/h) 

X lijsiX + (r + z)/2)l6/5|^(X + (z - r)/2)f/^\UX + {z + w)/2,r + w - z)|6/5 
1 

/l42/5 



1/5 



1/5 
1/5 



drdwdz ao{r / h)ao{w / h){V aQ){z / h) 

dX drdzdw ao{r / h)ao{w / h) {V ao) {z / h)\7ps {X + {r + z)/ 2) f 

dXdrdzdwao{r/h)ao{w/h){Vao){z/h)\7ps{X + {z-r)/2)f 

dXdrdzdw ao{r / h)ao{w / h)(y ao){z / h)\^s{X + {z + w)/2, r + w — z)\ 
< /^=^/'ll«o||^^||ao||f lll-aollf llV.lle'/'llellf < Gh^/' 
where we used the bound \\£,s\\2 ^ Gh from Proposition ll.2i The term Gg j can be bounded analo- 



-|3/5 
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gously. As for the term G^!., we write ^s{x, y) = is{{x + y)/2, x — y), and we notice that 
G^(X) = -^ f drdzdwaoir/h){Vao){z/h) 



h^ 



X ^,(X + (z - r)/2)e.(X + W + r/2, X + r/2)6(X + w + r/2, X + z- r/2) 
drdz ao(r//i)(yao)(z//i)Vs(X + {z - r)/2){l^i,){X + r/2, X + z - r/2) 



where ('^5Cs)(2^)2/) denotes the kernel of the operator £,^^3, i-e. 



{isis){x,y) = I dz^^{x,z)^siz,y) = / dz^six,z)^siz,y) 
(here we use the fact that, by definition, ^^ is symmetric, hke a^). Therefore 



\G' 



(l)||6/5 
3,s 116/5 



1 



drdzao{r/h){Vao){z/h)MX + {z - r)/2)(e,6)(^ + r/2,X - z - r/2) 

1/5 



6/5 



drdz ao{r/h){VaQ){z/h) 



X / (iXdrdzao(r//i)(Fao)(^//i)|V's(^ + {z - r)/2)f/^\{Us)iX + r/2,X + z - r/2)f/'' 



1 1/5 



< 



qqIIi II^q qIIi 

/l24/5 



1/5 



dX(irdzao(?-//i)(V^ao)(z//i)iV's(^ + (2 - r)/2)f 



2/5 



< 



dXdrdzao(?'//i)(V'ao)(2//i)I(^.6)(^ + r/2,X - z - r/2)|2 
|«o||f||ao||^^||yao||f||V>.||f(Tr(a.)2)3/5 /^ 

/l3/5 - '^'^ 



3/5 



where we used that 

— II*: Il4 <r Ch"^ 

from Proposition 11.21 The term G^!, can be controhed using the bound Tr (7^ — asOs) < Gh from 
Proposition 11.21 Since (7^ — a^Os) > 0, this implies immediately that Tr (7^ — a^ag)^ < Ch"^. 
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Therefore 



\G 



(l)||6/5 



< 



1 



4,s 116/5 - ^6 



1 

^1? 



dX / drdzaQ{r/h){Vao)(z/h){-is-ocsas){X + r/2,X + z-r/2) 

xV^,(X + (z-r)/2) 

nl/5 

drdz ao {r/h){VaQ){z/h) 

X / dXdrdz ao (r/h) {Vao)iz/h) | (7, - a,a,) (X + r/2, X + z - r/2) 1^/^ 



6/5 



1 1/5 



< 



«olli II^q qIIi 

/l24/5 



1/5 



x|V',(X + (z-r)/2)|6/5 

-|2/5 



(iXdr(izao(r//i)(yao)(z//i) lV's(X + (z - r)/2)|^ 



< 



f dXdrdz ao{r/h){Vao)iz/h) 1(7, - asas){X + r/2, X + z - r/2)p 
l«ollf ll«ofo^'ll^«ollf IIV'llf(Tr(7s-as«s)^)^/^ ^ ^,3/5 

/l3/5 - ^'^ 



3/5 



Finally, we compute 

11^(1) ||6/5 1 



5,s 116/5 - ^18/5 

1 
< 



dX 



/1I8/5 



drdz ao{r/h)V{z/h)js{X + r/2, X + z- r/2)i,{X + {z - r)/2, z) 
1/5 



6/5 



drdz ao(^/^)^(-2//i) 



X / dXdrdzao{r/h)Viz/h)\jsiX + r/2,X + z - r/2)f/^\UX + {z - r)/2, z)f/^ 

3/5 



< 



I l|l/5||T/||l/5 

IqqIIi \\v\\i 

/ll2/5 



dXdrdz ao{r/h)V{z/h)\-fs{X + r/2, X + z- r/2)\'^ 

2/5 



dXdrdz ao{r/h)V{z/h) |6 {X + {z - r)/2, z)f 

< /^^/^llaollf ||y||^^||ao||^^||V^||f (Tr7,^)^/^||6llf < Ch^' 
using that, by Proposition 11.21 Tr^'^ < h and ||^s||3 < ||6||hi < h^'^^'^. 

A Strichartz Estimates 

Let n > 3. A pair {r,q) € [l,oo) x [l,oo) is called admissible if 2 < r < 2n/{n — 2) and {2/q) 
n(l/2 - 1/r). Let 



D 



</./(*)= / dsU{t-s)f{s) 
Jo 

For any two admissible pairs (r, q) and {p, 6) we have 

I|0/IIl9([o,T],L'-(R")) < ^\\f\\L^'{[0,T],LP'{ 

where {5',p') are dual indices to {d,p) (i.e. 1/6 + 1/6' = 1/p + 1/p' = 1). In our analysis we use 
Strichartz estimates with (r, q) = (2, 00) and (p, 6) = (6, 2). In the latter case, the Strichartz estimate 
is referred to as the endpoint Strichartz estimate, and was proven in [iTj . 
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